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Abstract 

The magnetic properties of the molecular cluster Mni2Ac are due to the four Mn^+ ions which 
have spins S=3/2 and the eight Mn^+ ions with spins S=2. These spins are coupled by an exchange 
mechanism. We determine the four exchange couplings assuming a Heisenb erg- type interaction 
between the ions. We use exact diagonalization of the spin Hamiltonian by a Lanczos algorithm 
and we adjust the couplings to reproduce the magnetization curve of Mni2Ac. We also impose 
the constraint of reproducing a gap of 35 K between a S=10 ground state and a first excited state 
with S=9. We predict that there is an excited level with S=8 at 37 K above the ground state, 
only slightly above the S=9 excited state which lies at 35 K and the next excited state is a S=9 
multiplet at 67 K above the S=10 ground state. 

PACS numbers: 75.25.+z,74.10.Jin 



Electronic address: Nicolas . Regnault @lpmc . ens . fr 



^Electronic address: Thicrry.Jolicoeur@lpmc.ens.fi 



I. INTRODUCTION 



Molecular nanomagnets consist of a few paramagnetic ions coupled by exchange interac- 
tions and their properties are lying at the borderline between quantum and classical behavior. 
The best studied cluster so far is certainly Mni2Ac. At low fields and low temperature 
its ground state can be described as a single spin S=10. The reversal of this spin occurs 
via a macroscopic relaxation time at T=2 K. This discovery has prompted a detailed 
examination of its magnetic behavior because of the potential use for ultimate information 
storage. If we concentrate only on magnetic ions, then Mni2Ac is made of an external ring 
of eight Mn'^+ ions with spin S=2 and ring encloses a tetrahedron of four Mn^"*" ions with 
spin S=3/2. These spins are coupled by exchange so as to lead to a S=10 ground state. This 
means that there is a ferrimagnetic arrangement of the spins. It is also known that the S=10 
manifold is split by anisotropy into sublevels with — 10 < S*^ < +10. The properties of this 
subset of levels have been studied in detail by various experimental techniques because of 
the possibility of macroscopic quantum tunneling of the spin. The anisotropy responsible for 
the zero-field splitting of the S=10 manifold is smaller than the exchange interactions that 
determine first of all the ferrimagnetic ground state structure. For example the degenerate 

= ±10 states are separated from the states with = ±9 (belonging to the same S=10 
manifold) by a gap determined by inelastic neutron scattering ^ 14 K. So a reasonable 
strategy to study magnetic properties of this cluster is to ignore first any anisotropy terms in 
a model spin Hamiltonian by use of the simplest Heisenberg exchange and then, once these 
parameters are known, refine the treatment by including higher-order anisotropy terms. 

Recently the full magnetization curve M(H) of the cluster Mni2Ac has been obtained 0] 
by an experimental technique using explosive compression of the magnetic flux to access 
to very high fields in the megagauss range. Many if not all transitions between levels with 
different total spin have been measured. This kind of measurement is a direct probe of the 
inner magnetic structure of the cluster i.e. of the higher-energy scale couplings. 

In this paper we determine the values of the exchange couplings between the Mn ions by 
using a Heisenberg spin Hamiltonian to reproduce the high-field magnetization data. We 
also add as a constraint the value of the gap between S=10 and S=9. With these values we 
are able to predict the energies and multiplicities of some of the low-lying states above the 
S=10 ground state. This explains some features of the existing inelastic neutron scattering 
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data that were not understood previously. In section II, we expose the details of the model 
Hamiltonian and the numerical technique we have used. In section III, the details of the fit 
are given. Finally section IV contains our conclusions. 

II. THE SPIN HAMILTONIAN 

The core of the magnetic cluster of formula Mni20i2(CH3COO)i6(Il20)4 is depicted in 
Fig. There is an outer ring of eight ions Mn^"^ with spin S=2, S5, . . . S12 that encircles 
four ions Mn^^ with S=3/2, Si, . . . S4. In Fig. |^ we draw a simplified representation of the 
cluster taking into account only the connectivity of exchange interactions |0, |^ . There are 
four most important exchange paths (bonds in the following). They are Ji which relate 
each spin 3/2 to a spin S=2, there is also a J2 coupling between the outer ring and the 
inner tetrahedron involving four triangles i.e. eight bonds, then J3 between the four spins 
S=3/2 with a tetrahedral structure so there are six such bonds, and J4 between the S=2 
ions which means eight bonds. The simplest guess for the exchange is thus a Heisenberg 
spin Hamiltonian : 

H= Ji {Si ■ S7 + S2 ■ Sg + S3 ■ Sii + S4 ■ S5} 

+ J2 {Se • (Si + S4) + Sg ■ (Si + S2) + Sio ■ (S2 + S3) + S12 ■ (S3 + S4)} 
+ J3 {Si • S2 + S2 ■ S3 + S3 ■ S4 + S4 ■ Si + Si ■ S3 + S2 • S4} 

+ Ji {S5 ■ Se + Sq ■ S7 + S7 ■ Sg + Ss ■ Sg + Sg ■ Sio + Sio • Sii + Sn ■ S12 + S12 ■ S5](,I) 

where we use the labeling of the spins given in Fig. 2. In this work we ignore the effect 
of anisotropy and use the Hamiltonian Eq.(|I]). As a consequence the total spin is a good 
quantum number. The total dimension of the Hilbert space of the magnetic degrees of 
freedom is 10^ which is too huge for brute force diagonalization. Previous works [Q, ^ have 
used various approximate schemes usually based on some assumptions about the relative 
order of magnitude of the J's. In the Florentine coupling scheme [§] one starts from a large 
antiferromagnetic value for Ji. Hence the cluster is described as a first approximation by 
four dimers Mn^^-Mn^"*" with spin S=l/2 and four remaining Mn^"*" ions with spin S=2. It 
is then feasible to treat the system of four spins S=l/2 and four spins S=2. However it is 
not clear that this hypothesis is able to explain the experimental results at hand. There are 
many results from neutron scattering , magnetization ^ and heat capacity that point 
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to the presence of excited states not too far from the S=10 ground state. It is thus desirable 
to obtain values of the exchange parameters without recourse to any a priori assumption 
of their relative strength. A recent study by exact diagonalization [|13| using the point 
symmetry group of the magnetic cluster concluded that Ji=215 K, J2=85 K, J3=85 K, J4— 
64.5 K. Since these data refer to transitions involving the lowest energy levels for a given 
spin value, the Lanczos algorithm appears to be particularly well suited. 

The recent megagauss magnetization data @] provide further information on the strength 



of the exchange interactions. We have used the Lanczos algorithm |]10| applied to the 
Hamiltonian Eq.(|l]). If we use conservation of spin projection on a given axis, say it z, then 
the Hilbert space dimensions are much reduced. The largest subspace corresponds to S*^ = 
with dimension 8581300, while the subspace in which we have to find the ground state has 
dimension 817176 for 5*^ = 10. The corresponding dimensions are given in table |. Although 
the use of total spin would reduce further the dimensionality, it is much more complicated 
to program in an efficient way. 

In a given subspace we use the iterative Lanczos algorithm that brings the Hamiltonian 
into tridiagonal form. This is done typically in at most one hundred iterations which is 
enough to obtain the ground state energy with a precision of ~ 10^^°. The advantage of 
the Lanczos algorithm is that it requires only to perform the product of the Hamiltonian on 
a vector. Since spin Hamiltonians are very sparse matrices, one needs only to store in fact 
two vectors to use the Lanczos algorithm. Even with the huge dimensions that appear in 
table I, this is feasible on present day computers. 



If we now consider the effect of a magnetic field |Tl|, then due to the rotational symmetry 
of the Hamiltonian, the magnetic field couples to a conserved quantity and we just have to 
shift energies : 

Es^{B) = Es^{B = 0) - g^i^BS\ (2) 

where Es^{B = 0) is an eigenenergy of (0) in the sector with spin projection 5*^. As a 
function of the applied field there will be crossings of levels with different values of S^. 
The magnetization curve of the cluster Mni2Ac is a series of discrete jumps and plateaus 
due to the finite size of this spin system. It is only for the thermodynamic limit that one 
gets smooth magnetization curves for spin systems. The critical field corresponding to the 
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transition to 5*^ + 1 is given by : 

BcruiS' -^S' + l) = -^(Eo,s^+,{B = 0) - Eo,s^{B = 0)), (3) 

where Eq^s^{B = 0) is the ground state energy in the sector with given S^. We have used 
g=2 in this paper. 



III. REPRODUCING MAGNETIZATION DATA 

To determine the exchange parameters, we have used the data coming from various 
experimental sources. First it is well established that the ground state has total spin S=10. 
These states are exactly degenerate in our case since we make the simplifying assumption 
of exact rotational symmetry. Then the gap between this S=10 manifold and the first 
excited state with S=9 is known from magnetic susceptibility measurements [|1^ to be 35 K. 



The remaining piece of high-energy information comes from the megagauss experiment of 
ref. IQ. In this experiment one measures the differential susceptibility dM/dH vs. H. For 
a system with a series of discrete jumps as Mni2Ac, a spike in the differential susceptibility 
corresponds to a change of the spin of the ground state. The lowest-lying spike lies at 
Bi = 382 T and is interpreted ^ as the transition from the ground state with S=10 to 
an S=ll state and then there are three spikes at B2 = 416 T, B^ = 448 T, and B4 = 
475 T corresponding to the crossing of S=12, 13, and 14 states respectively. At higher fields 
there is a huge spike centered at 530 T which presumably corresponds to several unresolved 
crossings, maybe S=15 and S=16. Above this field value it is difficult to locate the other 
remaining transitions so no other values were determined in Ref. 0]. We have used the four 
values that are determined with good accuracy Bi, . . . , B4. We compute the corresponding 
theoretical values Bf^*^" and measure the quality of the fit by the following quantity : 

4 / R Rthco \ 2 

We have performed calculations of the levels with S^=9 up to S^=22 on a grid of values 
of the ratios of the exchange couplings J2/J1, J3/J1 and J4/J1 between +2 and -2 by steps 
of 0.1. In this range of parameters, we first reject values for which S^=10 is not the ground 
state. Then we look for regions with small e parameter defined in Eq.(^. In these regions, 
we fix the absolute scale of energy via Ji by requiring that the critical fields should be equal 



to the experimental value Bi = 382 T. This leads to approximate values for the other critical 
fields B2, . . . , B4. We next refine the search by including only the regions in which the gap 
S=10-S=9 is close to its experimental value of 35 K. This is now done by adjusting the four 
dimensionfull couplings Ji,...,J4. We observe that the biggest effect on the overall spectrum 
structure is due to Ji and J2. The best values are close to Ji 119 K, J2 ~ 118 K. This 
corresponds to e ~ 10~^. In Fig. ^ we plot the variation of the fitting parameter e in the 
J1-J2 plane close to the best values of these parameters. If we search for a region of small e 
by tuning these two parameters only then we find that this region is quite insensitive to the 
choice of J3 and J4. This is illustrated in Fig. ^ where e is given in the J3- J4 plane. 

Our next observation is that once Ji and J2 are determined, the remaining couplings J3 
and J4 may be varied to obtain more precisely the S=10-S=9 gap of 35 K. The agreement 
with the gap is given in Fig. |^ where we have plotted the gap value in the plane J3-J4. We 
find that J4 is significatively antiferromagnetic J4 23 K while it is difficult to give a precise 
estimate for J3 : it barely differs significatively from zero. A tentative value is J3 ~ —8 K. 
Taking into account the experimental uncertainties on the critical fields, we estimate that 
J3 and J4 are determined with an error bar of ~ 6 K. The determination of J3 and J4 does 
not affect much the values of Ji and J2, as can be seen in Fig. ^. 

As proposed in previous works, we find that Ji is larger than J3 and J4. However, 
we clearly need a second coupling J2 which should be close to Ji (taking into account 
the uncertainties on all the experimental results). The part of the spectrum relevant to 
magnetization data corresponding to this preferred set of parameters is given in Table II. 

It is interesting to compare our results with those already available in the literature. 
In the first works 0, the assumption of large Ji combined with a perturbative treatment 
lead to Ji=225 K, J2=90 K, J3=90 K, J4=0 K. In fact this set of parameters does not 
lead to a S=10 ground state, as first pointed out by Raghu et al. |T^. There is another 
set of parameters Ji=215 K, J2=85 K, J3=-85 K, J4=-45 K suggested by Chudnovsky [jl4|. 
However this set which has a correct S=10-S=9 ordering leads to a gap of 223 K, much 



too large. A recent study by exact diagonalization [T^ concluded that Ji=215 K, J2=85 K, 
J3=85 K, J4=-64.5 K. This set has a correct ordering of 9-10 levels with a gap which is 
adjusted to the experimental value of 35 K but we find that it does not lead to a satisfactory 
magnetization curve. This set gives critical fields equal to Bi = 192 T, B2 = 239 T, B^ = 
356 T and B4 = 406 T (the e parameter is 0.4 instead of our value of ~ 10~^). 
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With our preferred set of exchange parameters we can then compute some of the excited 
levels that are above the first excited manifold with S=9. We find that there is a S=8 
manifold sitting at 37 K above the S=10 ground state so very close to the S=9 states that 
have been used to constrain our fit. Then there is a S=9 multiplet which is found at 67 K 
above the ground state. This picture is close to what is found from susceptibility mea- 
surements where two close S=9 levels are necessary to reproduce the low-temperature 
behavior. However the scheme is not exactly the same : the multiplicities do not coincide. 

Our new picture has some interesting consequences for the interpretation of neutron data. 



There is some evidence from inelastic neutron scattering for a mode at 1.2 THz i.e. 70 K, 
this may correspond to the second S=9 multiplet that we find at 67 K. This is allowed by 
neutron scattering selection rule starting from the S=10 ground state. It is then normal 
that neutrons do not see the S=8 states at 37 K = 0.73 THz because of the selection rules 
however it not yet clear why neutrons do not see the first excited state with S=9 at 35 K = 
0.72 THz. 



IV. CONCLUSIONS 



We have studied the energy levels of a Heisenberg spin model Hamiltonian appropriate 
to describe the magnetic cluster Mni2Ac. We obtain a determination of four exchange 
couplings : Ji=119 K, J2=118 K, J3=-8 K, J4=23 K. Such a set of parameters reproduces 
the magnetization curve observed in megagauss experiments, leads to a S=10 ground state 
and a gap of 35 K to a first excited level with S=9 as measured experimentally. The 
numerical method we have used does not rely upon any approximations, so the main source 
of uncertainty in the values we quote comes from the measurement of the critical fields in 
the magnetization process. We predict that there is an excited level with S=8 at 37 K only 
slightly above the S=9 excited state which lies at 35 K and the next excited state is a S=9 
multiplet at 67 K above the S=10 ground state. These findings explain part of the existing 
neutron scattering data. 
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TABLE I: Hilbert space dimension of subspaces of given for the system of eight spin-2 and four 
spin-3/2. 
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TABLE II: Ground state energies in each fixed sector with corresponding critical magnetic fields 
for coupling constant set Ji=88.9 T 119 K), J2=88.0 T (« 118 K), J3=-6.0 T (« -8 K) and 
J4=17.0 T (ft! 23 K) for which the gap is 26 T (ft 35 K). 

Energy (T) Critical fields (T) 

10 -3196.2 

381.0 

11 -2815.2 



12 -2403.5 



13 -1959.4 



14 -1481.8 



15 -969.3 



16 -420.9 



17 164.0 



18 786.1 



19 1445.8 



20 2143.0 



21 2876.9 



22 3641.8 



411.7 



444.1 



477.6 



512.5 



548.4 



584.9 



622.1 



659.7 



697.2 



733.9 



764.9 
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Jl(T) 




FIG. 3: Variation of the fitting parameter e as a function of Ji and 32- The remaining couplings are 
taken to be J3=-6.0 T (-8 K), J4=17.0 T (23 K). The grey levels are proportional to - log(e x 10^). 
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FIG. 4: Variation of the fitting parameter e as a function of J3 and J4. The remaining couphngs 
are taken to be Ji=88.9 T (119 K), J2=88.0 T (118 K). The grey levels are now proportional to 
2/(1 + ex 10''). 
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FIG. 5: Variation of the gap as a function of J3 and J4. The remaining couphngs are taken to be 
Ji=88.9 T, J2=88.0 T. Grey levels are plotted via 1/1 + |(A - 26)/26|. 
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FIG. 6: Variation of the gap as a function of Ji and J2. The remaining couplings are taken to be 
J3=-6.0 T, J4=17.0 T. Grey levels are given by the relative variation |(A — 26)/26|. 
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